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1. Introduction 
Computer work, especially computer graphics, has been an important tool of 



^^ , Abstract. We describe a rigorous and efficient computer algorithm for build- 

. ■ ing a model of the dynamics of a polynomial diffeomorphism of C^ on its 

,S^ ' chain recurrent set, K, and for sorting points into approximate chain transi- 

tive components. Further, we give explicit estimates which quantify how well 
this algorithm approximates the chain recurrent set and distinguishes the chain 
transitive components. We also discuss our implementation for the family of 
Henon mappings, fa,c{^, v) = {^'^ + c — ay, x), into a computer program called 
Hypatia, and give several examples of running Hypatia on Henon mappings 

> 

(N 

O 

vQ ■ discovery in the field of dynamical systems. This paper is also concerned with the 

f^ I use of computers but it has a different goal. The goal of this paper is to rigorously 

^^ ■ establish some results on the location of the set of points where recurrent behavior 

^P I takes place. We focus on the class of polynomial diffeomorphisms of C^, which 

r^ ■ includes the widely studied family of Henon mappings. fa,c{x,y) = {x^ -\- c — ay , x) . 

We start by examining a rigorous and effective computer algorithm for building 
rj ■ a model of the dynamics of a map, /, on its chain recurrent set, 3?, and for sorting 

J~J I points into approximate chain transitive components. We call this algorithm the 

^ . box chain construction. The same basic algorithm has been studied previously, in 

'k^ \ different settings. Osipenko and Campbell ([20112]) approximate the chain recur- 

'V^ ■ rent set for a homcomorphism of a smooth, real, compact manifold. Eidenschink 

C^ ' ([5|) discusses a similar procedure for real flows. A philosophically related proce- 

dure is studied in [Hj |7] , but their case of interest is the attractor of a real map 
(rather than the chain recurrent set). |15| is a recent survey of this work. The focus 
in the previous body of work is to develop a very general procedure for rigorously 
approximating 3^ for continuous maps or flows in M". 

In this paper, we restrict our attention to polynomial diffeomorphisms of C^, 
which allows us to adapt the box chain construction and its implementation to be 
more efficient. In addition, we establish estimates on the accuracy of our model. 
These are explicit estimates, involving only the inputs to the program, which quan- 
tify how closely our box chain construction approximates the chain recurrent set. 
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This allows us to predict when an execution of the program will successfully sep- 
arate the distinct chain transitive components. We contrast this with the work of 
Dellnitz and Hohmann ([HI), which gives results on the accuracy of the approxi- 
mation in the case that the map is hyperbolic; however, their estimates involve 
constants of hyperbolicity, which they do not discuss how to calculate. Osipenko 
and Campbell (|22) S^^^ approximation estimates, in terms of some of the inputs to 
the program as well as selected output, hence the accuracy of their approximation 
can only be measured after execution of the algorithm. 

Our first result on accuracy of the box chain construction for polynomial diffeo- 
morphisms of C^ can be paraphrased as follows. 

Theorem 1.1. Let f be a polynomial diffeomorphism 0/ C^ (of dynamical degree 
d > 1). Suppose "Bq is a closed box in C^ containing the S^-chain recurrent set, 
51{6q), for some S'q > (for example, take 23o, and 6 ~ S'q > as in Provosition \2. b]) . 
The box chain construction produces sequences of constants, {5n\ and {£«}; di- 
rected graphs, {r„}, and regions in C^, {^n}, for n>\, such that 

(1) (5„ ^ e„ and both e„ j and (5„ i as n ^ 00, 

(2) the vertex set ofT^ is a collection of closed boxes {i?^} in C^, which have 
side length at most En, 

(3) each !B„ is the region in C^ defined by the union of the B 

(4) CBn} is nested, i.e., {• • • C Sn C • • • C Si C Sq}, 

(5) there is guaranteed to be an edge in r„ from B^ to B" if f{B^) intersects 
a Sn-neighborhood of B^-" , and 

(6) we can calculate explicit sequences {e^}^i and {(5,'j}5^i, and an explicit 
constant C , such that for every n > 1, 

(a) e^ I 0, in particular e„ < e^ < (5„ + Ce„, 

(b) (5^ is nonincreasing and converges to zero, with 5^ < 5n, and 

(c) ^8'^) C S„ c 3?«). 

Definition 1.2. For any n>\, suppose (e, 5, F, 53) = (e„, (5„, F„, B„) are produced 
by the box chain construction at step n, and satisfy Thcorcm ll.il 

We call the region S an (e, S)-box chain recurrent set, and the graph F an (e, S)- 
box chain recurrent model of f. Each edge-connected component F' of F is called 
an {e,5)-box chain transitive component. 

Conclusions (1) through (5) of Thcorcm ll.ll foUow immediately from the descrip- 
tion of the box chain construction, given in Section 12.21 Conclusion (6) is our first 
significant a priori estimate on the accuracy of our model, and is established in 
Sections 12.31 through 12 . 61 There we show that for the case of Henon mappings, for 



k' 



any chosen R' > 1 -f |a| -f ^(1 + \a\)'^ + 4 |c| /2, we get 

S', = [iRr-il + \a\)R'-\c\]/2, 

5[, = min(<5;, [-(2i?' + \a\ + 1) + ^J{2R' + |a| + 1)2 + 4,5„] ) 

4 = 5n+en{l + \a\ + 2R' 



, and 



_2 



A box chain recurrent model of / satisfies the definition of a symbolic image of 
/, given by Osipenko in ^^. Osipenko and Campbell (|U) derive similar estimates 
to (6), but their version of ej^ depends on measuring the size of the images of boxes 
Bfe, after they have been computed. 
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We have implemented our efficient box chain construction for Henon mappings 
into a computer program we call Hypatia. Using Hypatia we have produced many 
examples of box chain recurrent models of Henon mappings. ^ The box chain 
construction also has an immediate analog for polynomial maps of C, which we 
include in Hypatia for quadratic and cubic polynomials. We keep the arithmetical 
computations in Hypatia rigorous using interval arithmetic with directed rounding. 
This method was recommended to us by Warwick Tucker, who used it in his recent 
computer proof that the Lorenz differential equation has the conjectured geometry 
f|24|'). See Appendix IbI for a brief introduction to interval arithmetic. 

The examples we are most interested in studying with Hypatia are Henon map- 
pings which are hyperbolic. These are the class of maps which have the "simplest" 
chaotic dynamics, and are in fact stable under small perturbation. Thus hyperbolic 
maps are the most amenable to rigorous computer investigation. In fact, Bedford 
and Smillie ([H|) have shown that for hyperbolic polynomial diffeoniorphisms of C^, 
the chain recurrent set is well-behaved, in that it consists of simply the Julia set 
together with finitely many attracting or repelling periodic points. 

The simplest hyperbolic Henon mappings can be described in terms of the dy- 
namics of some quadratic polynomial. In fact, if fa,c is a Henon mapping with a 
sufficiently small and c is such that the polynomial Pdz) = z^ + c is hyperbolic, 
then /I J is topologically conjugate to the function induced by P on the inverse 
limit lim^(J, P) f|14|). In this case, we say that / is described by P, or simply that 
/ exhibits one dimensional behavior. The work of Hubbard and Papadontonakis 
(^IP), and more specifically the work of Oliva (CHl)i suggests points in parameter 
space which are conjectured to be hyperbolic and to have interesting properties. At 
the moment though neither the hyperbolicity nor the interesting properties have 
been established rigorously. 

Thus a significant problem in the study of the Henon family is to understand 
which maps exhibit one dimensional behavior, and to describe the behavior of 
maps which do not. Motivated by this problem, in this paper we use the box chain 
construction and its implementation in Hypatia to build box chain recurrent models 
for several interesting examples of Henon mappings. Further, we use the results of 
this paper as the first step in a study of the property of hyperbolicity for polynomial 
diffeomorphisms of C^ in |11| . and in a study of hyperbolicity {i.e., expansion) for 
polynomial maps of C in |12j . |1()| contains an earlier version of the work of this 
paper, as well as that of |11II12J . 

Example 1.3. One of the simplest Henon mappings which appears not to exhibit 
one dimensional behavior is fa.c with (a, c) = (.125, —1.24). Oliva (jEl) gave combi- 
natorial evidence that this diffeomorphism is hyperbolic with a period two attract- 
ing cycle, but is not described by a quadratic polynomial. Using our program Hypa- 
tia, we applied the box chain construction to the nearby map (a, c) = (.15, —1.1875), 
which seems to be topologically conjugate to {a,c) — (.125,-1.24). We computed 
a sequence terminating at a box chain recurrent set Sg. 

For a qualitative estimation of the accuracy of a box chain recurrent set, ®, 
we can sketch the intersection of "B with a dynamically significant one-dimensional 
submanifold of C^: the unstable manifold of a saddle fixed point, which has a 
natural parameterization by C. This process is explained in Appendix I A. 31 and 
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Figure 1. For fa,c the Henon mapping with c = -1.1875, a ~ .15, 
m appears to be J and a period 2 sink. Shown in this figure are two 
box chain recurrent sets, restricted to the unstable manifold of a 
saddle fixed point, with its natural parameterization. On the top, 
boxes are of side length 2R/2^ and 2R/2'^, where R = 1.9. On the 
bottom is a refinement obtained from subdividing once some of the 
boxes on the top. The innermost, darkest regions in these figures 
are contained in the box chain transitive component containing 
the 2-cycle. Each box chain transitive component containing ,/ is 
shaded two tones by a heuristic algorithm, in order to illustrate 
how close the component is to J. The top is the crudest box 
chain recurrent model which separates J from the sink. In the 
bottom figure, the small dark spots skirting the inner circles of the 
neighborhood of J form a box chain transitive component of points 
which arc pseudorecurrent but not recurrent. 



Section 15^ In Figure ^ we use this procedure to illustrate two box chain recurrent 
sets, "Bt and Sg, for the Henon mapping with {a,c) = (.15,-1.1875). 

Assuming the conjectural dynamics holds, the box chain recurrent models con- 
structed for Example 1 1.31 are both successful in the sense of the following. 
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Figure 2. A box chain recurrent set for the Hcnon mapping /^ c, 
with a = .3, c = —1.17, restricted to the unstable manifold of a sad- 
dle fixed point, with its natural parameterization. Here boxes are 
of side length between 2R/2'^ and 2i?/2^, where R = 2.01. Lighter 
gray points were heuristically found to be in A'+. Unfortunately, 
this box chain recurrent model does not separate J from either the 
fixed sink or the attracting three-cycle. 



Definition 1.4. Let / be a polynomial diffeomorphism of C^ (of dynamical degree 
d > 1). Let r be a box chain recurrent model of /. We call F separating if there 
are two chain transitive components, Jl^ and Jl^ (of 3^), which lie in different box 
chain transitive components of F. In this case we say F separates Jl^ and Jl'^. 

Further, we call F fully separating if it separates every pair of chain transitive 
components. 

Example 1.5. Another interesting example studied by Oliva (jIH]) is the Hcnon 
mapping fa.c with (a, c) = (.3, —1.17). We call this the 3-1-map, because it appears 
to be hyperbolic, with 31 consisting of three chain transitive components: J, an 
attracting fixed point, and an attracting cycle of period three. In contrast, quadratic 
polynomials cannot have more than one attracting cycle, thus this map appears 
not to exhibit one dimensional behavior. We applied the box chain construction 
to the 3-1-map, but were surprised to be unable to obtain a separating box chain 
recurrent model. The best box chain recurrent set we obtained is shown in FigureEl 
intersected with the parameterized unstable manifold of a saddle fixed point. 

Our difficulties with the 3-1-map motivated the following theorem, in which we 
calculate explicit bounds on e and 6 to guarantee that an (e, 5)-box chain recurrent 
model will separate the fixed sink from the 3-cycle and the Julia set. This gives a 
theoretical quantification of the computational difficulty of studying the 3-1-map. 

Theorem 1.6. Suppose fa^c *s a Henon mapping with an attracting fixed point p, 
with Ai 7^ A2 eigenvalues of Dpf, and A = niax(|Ai| , IA2I). Set 



IA1-A2I' 



(2+|Ai| + |A2|)(2 + A2 + |a|)- 



Let F be an (e, 5)-box chain recurrent model of f . Let M > 1 satisfy 6 < e/M . 
Set 

K=[l + l/AI + max{l, (1 - A) V7 + 2 ||pi| + \a\}] . 
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If e < ^ ( ^'* + \/k^ + t(1 — A)2 ) , then T separates the fixed sink from every 
other chain transitive component of 51. 

This theorem apphed to the 3-1-map yields that boxes of side length less than 
4 X 10~^ would guarantee separation. However this is several orders of magni- 
tude smaller than current resources allowed us to achieve with Hypatia. This 
demonstrates the need for the development of a more sophisticated construction 
for rigorously approximating chain recurrent sets of complex Henon mappings. 

Finally, we outline the remaining sections. In Section [3 we describe the box 
chain construction for polynomial diffeomorphisms of C^ , and prove Theorem 11.11 
by calculating explicit estimates on how well a box chain recurrent set approximates 
the chain recurrent set. In Section Owe use some dynamical information about the 
map to develop two enhancements to the basic construction, significantly improv- 
ing computational efficiency. In Section 0] we show a theoretical limitation of the 
box chain construction by proving Theorem 1 1.61 and applying the estimates of the 
theorem to the 3-1-map. In Section[Slwe discuss examples generated with Hypatia, 
for Henon mappings and a polynomial map of C We have included relevent back- 
ground material on the chain recurrent set and the dynamics of Henon mappings 
in Appendix fXl In Appendix IBI we sketch the basics of interval arithmetic. 

Acknowledgements. We thank John Smillie for providing guidance on this 
project, John Hubbard, Greg Buzzard, and Warwick Tucker for many helpful con- 
versations on the topic, James Yorke, John Milnor, and Eric Bedford for advice 
on the preparation of this paper, Robert Terrell for technical support, and Michael 
Benedicks for pointing out to us that |15j describes a procedure similar to ours. 

2. The box chain construction for polynomial diffeomorphisms of C^ 

In this section we start by giving an outline of the box chain construction, then 
discuss how we carry it out for polynomial diffeomorphisms of C^ , to calculate the 
estimates leading to Theorem ll.il To calculate our estimates, we assume a polyno- 
mial diffeomorphism of C^ (of dynamical degree d{f) > 1) is a finite composition of 
generalized Henon mappings, which are maps of the form f{x,y) = (p(x) — ay,x), 
for p nionic of degree greater than one (see Appendix lA)l . 

2.1. Efficient neighborhoods. Before we begin our theoretical calculations, we 
want to specify that we do not use the euclidean metric. It is more natural for 
computer calculations to consider vectors in M^", rather than C", and use the L°° 
metric, rather than euclidean. Thus throughout this paper, ||-|| will denote the L°° 
norm on R^", so that for a vector x = {xi, . . . , Xn) & C", 

||a;|| = max{|Re(xfc)| , |Im(a;fe)| : 1 < fc < n}. 

Also, let y^lS, r) denote the open r-neighborhood about the set S in the metric doo 
induced by \\-\\, e.g., 

?^(0,r) = {x=(xi,...,a;„)eC": |Re(2;fe)| < r and \lm{xk)\ < r}. 

We use the simpler notation |-| for dimension n = 1. This metric is uniformly 
equivalent to the euchdean metric on C",||-||g, since -^ ||a:;||g < ||x|| < |ja;||g. 
Neighborhoods are slightly different with respect to two uniformly equivalent norms, 
but the topology generated by them is exactly the same, thus they can practically be 
used interchangeably. Similarly, the e-chain recurrent set 3i(e) depends on choice of 
metric, but since 3?'^(e) C 3?(e) C 'H'' {^/2ne) , the chain recurrent set 31 = ne>o3?(e) 
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is the same for any metric uniformly equivalent to euclidean. Thus throughout we 
use 3^(e) as defined by our norm. 

Remark 2.1. When we say 60a;, we mean a ball around a point in this norm. Note 
a box is also a vector of intervals, so boxes are neighborhoods which are easily 
manipulated with interval arithmetic. 

2.2. The box chain construction. As suggested by the statement of Theo- 
rem ^^ the box chain construction is an inductive process. We use the idea that 
any 31(e) consists of precisely the e-pseudopcriodic orbits. Below is an outline of 
the basic construction. 

We start with a polynomial diffcomorphism / of C^ (of dynamical degree d{f) > 
1), and let F be any interval extension of /, i.e., i^ is a function on interval vectors 
such that for any box B, the image F{B) is a box containing f{B) (see Appendix IbI 
for background on interval arithmetic). 

(0) Given / and an interval extension F, choose a small constant 6'q > 0, and a 

closed box So in C^ which contains ^{5q) (of /). Let eo be the side length 

of the box Sq. Choose Sq such that eo ^ <5o > ^o- 
(n) Let n > 0. Suppose "Bn is a closed region in C^, consisting of a collection 

of boxes V„ = {-Bfe } of side length at most e„, and such that Bn D 3i((5^), 

for some S'^ > 0. Suppose (5„ > S'^ is given (if n > 1, then (5„ is given by 

step (?7. — l)-(ii).) 
(i) Equally subdivide the boxes in V,i. That is, choose m > 1, set e„+i = 
En/iTi, and place a grid of m'^ boxes inside each box of V„. This defines 
a new collection, W„+i, in which each box has side length e„+i. 

(ii) Build a graph approximating the map / on Wn+i. Specifically, choose 
some 6,1+1 such that Sn+i < 6n/2 and Sn+i <C Sn+i (for example, 
for every n > 0, set (5„ = e„/1000). Then compute a directed graph 
T„+i whose vertices are the boxes in Wn+i, and such that there is 
guaranteed to be an edge from box -BJ^^^ to box 5"+^ if F{B'^^^) 
intersects a 5„+i-ncighborhood of B"^^ . 

(iii) Find the subgraph of T„_(.i consisting precisely of the vertices and 
edges which lie in cycles. Call this subgraph r„-|_i . Let Vn+i = {B^"''^} 
be the vertices of r„+i, and define 'Bn+i as the union of the boxes in 
V„+i. 

Remark 2.2. The box chain construction immediately implies that statements (1) 
through (5) of Theorem II. II arc satisfied. 

Remark 2.3. The only difference between this basic box chain construction and 
the procedure used by Osipenko (and others as discussed in the introduction) is the 
presence of the constants S'^ and Sn- In order to approximating 31, these constants 
are uneeded, and can all be taken to be zero. However, in order to apply this 
construction to the problem of proving hyperbolicity, as we do in ^1^2i positive 
(5's are essential. 

The first step in verifying the usefulness of this construction is to show that 
if 33„_i D Jl{S'^_i), then there exists a S'^ such that !B„ D '^{^'n)- This fact for 
polynomial diffeomorphisms of C^ follows from Theorem 12.151 and Theorem 11.11 
proved at the end of this section. 
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2.3. Trapping Regions for 51. In order to apply the box chain construction to 
polynomial diffeomorphisms of C'^, we must first provide the base case, i.e., step 
(0) above. For this we prove Proposition 12.61 in which we calculate an explicit 
trapping region !Bo for the 5-chain recurrent set of polynomial diffeomophisms of 
C^. In particular, given a map / and ^ > 0, we give explicit R' such that 3?((5) C 
'Bo^{\x\<R',\y\<R'}. 

First we quantify how, for polynomial diffeomorphisms of C^, infinity in the 
X direction is attracting, while infinity in the y direction is repelling for /, and 
vice- versa for /^^. A version of the following lemma is given in |5] (also see \2'6\). 

Lemma 2.4. Let f be a polynomial diffeomorphism of C^, with d(f) > 1. Let 
(5 > 0. Then there is an R > 1 and an R' > R, such that if \x\ > R' and \x\ > \y\, 
then\\fix,y)\\>\x\+2S. 

If fa c is a Henon mapping, then R = i(l + |a| + ■\/(l + \a\Y + 4 |c|) and 25 ~ 
{Rr-{l + \a\)R'~\c\. 

Proof. Assume / is a generalized Henon mapping, f{x, y) = (z, x) = {p(x) — ay, x), 

withd =deg(p) > 1. Sop(x) = x'^+Cd-icc'^^^H I-cq. Let q(r) = r''--\cd-i\r'^~''^- 

■ ■ ■ — I Co I — (1 + laD^"- Then there is an i? > such that q is monotone increasing 
on [R, oo), with q{R) = 0. Note if r > 0, then q{r) <r'^ ~r, thus i? > 1. 

Since q \s a polynomial and is monotone increasing on [i?, oo), with q{R) ~ 0, 
we see q : {R, cxd) -^ (0, oo) is invertiblc. Thus given 6, define R' £ {R, oo) by 
q{R') = 25. 

Let {x,y) e C satisfy |.t| > R' and |.t| > \y\. Then 

\z\ > \p{x)\ - \a\ \y\ > \pix)\ - \a\ \x\ > q(|x|) + \x\ > 25 + \x\ . 

If / = fm°' ■ 'o/i, then each composition moves \z\ farther away from |x| additivcly. 
Thus, let R'^, satisfy qk{R'j^) = 25/m and take R' = max{i?fc}, for 1 < k < m. 

Note if fa.c is a quadratic Henon mapping, fa,c{x,y) = (x^ + c — ay,x), then 
q{r) = r^ — \c\ + (1 + \a\)r, hence we easily achieve the claimed bound. D 

Definition 2.5. Let 5 > 0. Following [2], we define the "trapping regions", for R' 
given by Lemma 12.41 by: 

Bo = {|a-|<i?'and \y\ < R'}; 
S„ = {|a:|>i?'and \x\ > \y\}; 
23+ = {\y\>R' and \y\>\x\}. 

Note the regions "Bq — 'Bo{5) depend on 5, but in each instance the 5 will be clear 
from context, so we will rarely use the notation 'Bq{5), rather we simply use 530. 

The sets V,V^ introduced in [2] are equal to !Bo,23q for (5 = 0, and satisfy 
K+ CVUV+, K~ CVUV-, and A' C y (0). Choosing R' larger than R 
allows us to preserve these relationships and trap (5-pseudo-orbits as well. 

Proposition 2.6. Let f be a polynomial diffeomorphism ofC'^, with d{f) > 1. Let 
(5 > and let 23o be as in Detinition \2.5\ Then 3^((5) C So. 

Proof. Given Lemma Ol we have /CBg ) C ®^ and /CBg ) n N(So, 2d^) = 0. Thus 
if p G Sq , then p is not in 3^(5), since the images f{xk) move by at least 25 in the x 
direction, and so the Xk+i coming back in by only 5 makes it impossible for a;„ = p. 
Similarly, for p S 53^ look at the chain backwards to contradict J-chain recur- 
rence. D 
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To get an idea of the size of R, note that R = 2 for c = 2. Since the Mandelbrot 
set is contained in {x : |a;| < 2}, for the parameters we tend to study we have 
1 < i? < 2. For a Hcnon mapping, the values of R are also close to this range. 

2.4. Defining tiie grapiis T and F. In steps (ii) and (iii) of the box chain con- 
struction, we compute graphs T and F representing the action of / (or F) on our 
collection of boxes. The following terminology will ease our discussion of these 
graphs. 

About notation: if F is a graph, then V(F) denotes the vertex set of F, and 
£(F) denotes its edge set. Also wc often discuss one subset of a collection of boxes 
V = {Bk}^^i at a time, and since the ordering is unimportant we avoid double 
subscripts and simply use {Bq, _Bi, . . . , Bi}. 

Definition 2.7. Let T = T„ be the directed graph built in step (n)-(ii) of a box 
chain construction. Then we say T is an (e,d)-box chain model of F. 

In addition, when we say T is an (e,(5)-box chain model of /, we mean the 
"theoretically ideal" model, using the ideal interval extension of /, i.e., require 
F{Bk) = Hull(/(Sfe)), for aU boxes Bfc. 

Note for any interval extension F of /, wc know F{B) D Hull(/(i?)) for any box 
B. Thus any result which is true for all interval extensions _F of / is also true for 
/. Thus our default is to discuss box chain models of F, and only use box chain 
models of / when trying to be precise about theoretical estimates. 

Suppose T is an (e, (5)-box chain model of F. Note by DefLnition ll.2l the subgraph 
F consisting of the cycles of T is called an (e, (5)-box chain recurrent model of F. 

Remark 2.8. Both a box chain model and a box chain recurrent model of / satisfy 
the definition of a symbolic im,age of f, given by Osipcnko in |19| . 

When the context is clear, or the distinction is unimportant, we simply refer to 
a box chain recurrent model as a box chain model. Thus we use the symbol T for 
any graph which is either a box chain model or a box chain recurrent model, and 
reserve F for box chain recurrent models. 

The following standard concept in graph theory will help us to analyze F. 

Definition 2.9. Let T be a directed graph. If there is a path from vertex v to 
vertex u, then we say u is reachable from v. A strongly connected com,ponent (SCC) 
is an equivalence class under the "are mutually reachable" equivalence relation. If 
an SCC consists of only one vertex, it must have an edge to itself. 

Note the following easy relationship. 

Lemma 2.10. Let T he any directed graph. Let F he the subgraph of T consisting 
of the vertices and edges which lie in cycles. Then F ~ (F^UF^U. . .UF') is precisely 
the subgraph of T consisting of the union of the SCO's of T. 

Recall from Definition lI.2l that the box chain transitive components of F are de- 
fined as the edge-connected components of F, hence these are precisely the SCC's of 
T. The relationship between a box chain transitive component and an e-chain tran- 
sitive component is made precise in Corollarv l2.17l and justifies our terminology. 
Thus the decomposition of a graph F into its SCC's is analogous to the partitioning 
of the chain recurrent set into it's invariant pieces, the chain transitive components. 

gives a standard procedure for decomposing a graph into its SCC's. 
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2.5. Estimates on T and F. Now we can calculate estimates which quantify our 
models. Wc begin with a lemma on the size of the image of the boxes. 

Lemma 2.11. Let T„ be an (sn,Sn)-box chain model of f. Then there exists 
r„ > (depending on e„, /, and "Bo) such that for e'^^ = (r„ + l)£,i + Sn, and for 
any Bk G V(T,i), we have: 

(1) the side length of the box Hull{f{B}S)) is less than or equal to Snrn, and 

(2) if {k,j) e £(T„), then for any Xk G Bk and any Xj € Bj, \\f{xk) - Xj\\ < 

For fa^c o, Henon mapping, we may take Tn = e,i + (2i?' + \a\), where R' is as in 
Proposition \2.6\ 

Proof. The second item follows immediately from the first item, and the fact that 
if Bj nJ^(Hu\l{f {Bk}) , S) ^ 0, then there must be an edge from Bk to Bj. We prove 
the first item using the linearization of / to approximate it. 

We assume / is a generalized Henon mapping, so f{x,y) — {pix) — ay,x), p 
monic of degree d > 1. Let B G V(T) and {z,w) G B. The linearization of / at 
(z, w) is 

Lzf{x, y) = f{z, w) + DJ{x - z,y~w) ^ {p{z) - p'{z){x - z) - ay, x). 

Note Lzf{z, w) = f{z, w), for any w. Next, we observe 

\\f{x,y)-LJ{x,y)\\^\Mx)-p{z)-p'iz){x-z),0)\\<Y,\-^^^\x-zf. 

fc=2 

If {x,y) is also in B, then |a; — z| < e„. Since p is a polynomial, and z G So = 
J^(0, i?'), then for any fc > 0, there exist Tk > such that \p^''\z)\/k\ < Tk. Hence, 

k=2 ■ fc=2 

Now we need to bound 

\\Lzf{x,y) - L^f{z,w)\\ = nmx{\x~ z\,\p'{z)(x- z)-a{y-w)\). 
If {x, y) is in _B, then wc also know \x — z\ < Sn, and 

\p'{z){x - z) - a{y - w)\ < Tie„ + \a\ e„. 
Finally, we put the above pieces together to compute, for any (a;, y), (z, w) in B, 

d 

11/(2-, y)-/(z,uO||<^TfeeJ^ + £„ max (l,ri + |a|). 

k=2 

Hence, if r„ is set to be the maximum above divided by e„, then we have that 
the diameter of the set f{Bk) in the metric ||-j| is less than or equal to r„e„, hence 
the side length of Hull(/(i?fc)) is less than or equal to r„£„. 

For fk a generalized Henon mapping, we have rn{fk) as defined above. Now, 
if ./ = .fm°---° /i, then set r„ = r„(/i) • • •r„(/m). This suffices, for if we let 
s{fk,£) = rn{fk)£ be the bound on the diameter of a box of size e under fk, then 
we get si = s(/i,e„) = e„r„(/i), then S2 = s(/2 o/i,e„) = s(/2,si) = sir„(/2) = 
£rnUi)Tn{f2), ctc, and s„ = s(/„ o • • • o /i) = s(/,„, s™_i) = er„(/i) • • • r^ifm)- 

In the case of /q ^ a Henon mapping, since p{z) ~ z"^ + c we can compute that 
Ti = 2R' and T2 = 1, and recall from Proposition 12.61 we have it!' > 1. Hence 
rn = en + max [2R' + \a\ , 1) = £„ + (2i?' + \a\). D 
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This lemma give us the foUowing component of Theorcm ll.il 

Corollary 2.12. Item (6)-(a) of Theorem M.lM s satisifed. 

Proof. First, we show there exists a C > such that for any ri > 1, r„ < C — 1, 
hence e'^ < Sn + Csn ■ 

Suppose for all n > 1 that £„ < 1 . Then we know ef^ < e„ for any fc > 1 . Define 
C = I + X]a:=2 -^/^ + max(l,Ti + \a\) for a generalized Hcnon mapping, or in the 
manner analogous to the previous proof for a composition {i.e., for / = fm°- ■ -o/i, 
take C = (C,„ - 1) • • • (Ci - 1) + 1). Then we see r„ < C - 1. 

On the other hand, since e„ is a decreasing sequence (see Remark |2 .21) . at most 

En > 1 for < n < iV. Thus if we set C = 1 + ma.x{ri, . . .rN,J2t=2Tk + 
max(l,ri + \a\)}, then r„ < C - 1 for all n > 1. 

Finally, the above lemma defines e',^ so that e^ > e„. Hence we have e„ < e^ < 
5n + Csn- By construction (see Remark |2.2f) we know both £„ and (5„ decrease to 
zero as n — > oo, hence e'j must as well. D 

The next lemma will allow us to find 5' for which a box chain recurrent model 
traps all 5' pseudo-periodic orbits, i.e., 3l((5'). 

Lemma 2.13. Let T„ be an (sn,Sn)-box chain model of F. There exists an ij = 
ri{n) € (0,(S„), such that if x g J4{Bk,ri) and f{x) £ J4{Bj,ri), then there is an edge 
from Bk to Bj in T„, i.e., {k,j) £ £(T„). 
For fa.c a Henon mapping, we can take 

1 



il{n) = - [-i2R' + \a\ + 1) + ^i2R' + \a\ + If + ASr. 



where R' is as in Proposition 

Proof. Again we give the proof for f{x, y) ~ ip{x) — ay, x), p monic of degree d > 1. 
Let {x,y) G 'N{Bk,ri) and f{x,y) G 1^{Bj,t]). Let {z,w) be a point in Bk which 
realizes this minimum distance (since Bk is closed), i.e., \\ {x, y) — (z, w)|| < rj. Then 
in order to guarantee that (fc,j) S £(T), we just need \\f{x,y) — f{z,w)\\ < (5„ — rj- 
But examining the proof of Lemma 12.111 since (z, w) G Bk, we see that we have 

d 

\\fix,y)-f{z,w)\\ <^Tfer;^'+7?max(l,ri + |a|). 

fc=2 

Thus, we just need rj to satisfy 

d 

J2 '^^ri^ + ^ (^^^ (1' ^1 + 1^1) + 1) ~ '^" < 0- 

k=2 

Let g„(t) = J2k=2'^kt'' + tmax(2,Ti + \a\ + 1) - 5,,. Then we set rj = 77(71) to 
be the smallest positive root of q„ . If / = /,„ ° ■ ■ ■ ° fi, we may take rj = 
min{?/(/i), . . .,Tj{fm)}- For fa.c a Hcnon mapping, we get 9„(t) ^ t^ + t{2R' + 
\a\ + 1) — Sn, which leads to the claimed bound. D 

Corollary 2.14. Recall "Bq denotes the initial bounding box, with ^{6'q) C Sq for 
some S'q > 0. For each n > 1, set S'^ = min(77(7i), S'g), with rj(ri) as in Lemma \2.1!A 
Then item (6)-(b) of Theorem \l.l\ is satisfied. 
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Proof. Lemma 12.131 states that 77(71) < (5„, hence S'^ < Sn- Next, note from the 
proof of Lemma 12.131 that 77 is a monotone increasing function of i5„ > 0. By 
R,emark l2.2l thc sequence Sn decreases as n increases. Hence as n increases, we get 
r]{n) decreases. Thus (5^ is at most equal to S'q for a finite number of terms N, and 
then for n > iV, (5,'^ decreases. D 

2.6. Quantifying the accuracy of an (e, (5)-box chain recurrent set S. In this 

section, we combine the lemmas of the previous section to achieve Theorem 12.151 
which immediately implies Theorem 11.11 along with several relevant corollaries. 

Theorem 2.15. Suppose T„ is an (£„,5„)-6oa; chain model of f, and r„ is the 

iSn-,Sn)-box chain recurrent model of f consisting of the SCO's o/T„. Let e[-^ he as 
in Lemma \2. 1 1\ and let S'^ be as in Corollary \2.14\ 

Let 23 (r„) = "Bn denote the region in C^ covered by the box vertices ofT^- Then 
item (6)-(c) of Theorem \l.l\ is satisfied, i.e., 

no c s(r„) c ne'n)- 

Proof. We use the fact that precisely the vertices which lie in cycles in T„ are those 
which are in some strongly connected component of T„. 

We handle the two inclusions separately. First consider the inclusion 3?((5,^) C 
"BiTn)- We first need to establish 3?((5^) C ^(Tn), which we prove by induction. For 
n = 0, we can choose any 5'^ > 0, and produce the box Sq such that 3?((5q) C Sq by 
Proposition 12. 61 Consider Tq = Fq to be the graph with the single vertex 'Bq with 
an edge to itself. Now suppose n > 1, and assume Ji{S'^_i) C S(F„_i). Observe 
that since the boxes of T„ are obtained from subdividing the boxes of F„_i, we have 
'B(T„) = •B(F„_i), hence 3?(<5;,_i) C ■B{Tn). Also, since 5'„ < j;_i (CoroUaryEH 
we have 5i{S'„) C 3^(^;_i). Thus 5i{S'„) c S(T„). 

Suppose p G 3?((5'„). Then there exist xi ~ P,X2, . ■ . ,x,„_i,x,„ = p such that 
\\f{xk) — Xk+i\\ < S'^ for 1 < fc < TO. Note that Xk € 3^(<5^), for 1 < A: < to. Hence 
each Xfc e 23 (T„) as well. Then there are boxes Bk € T„ such that Xk G Bk for 
1 < fc < TO. Since \\f{xk) - Xk+i\\ < 6',^, we have f{xk) e 'N{Bk+i,S',J. Since 
S'n ^ ^n, there is an edge in T„ from Bk to Bk+i- 

Hence, p is in a box Bi which lies in a cycle of T„, Bi —t B2 —* ■ ■ ■ ^ B^-i -^ 
Bi. Thus Si e V(F„), hence p € Bn- 

For the second inclusion, 23 (F„) C 3^(eJJ, suppose p G 23 (F„). Thus p lies in some 
box Bi which lies in a cycle Bi ^ B2 ~^ ■ ■ ■ ~* Bm-i — *■ Bi in T„. Recall e„ is the 
side length of the boxes in the graph F„. Let xi = x„i = P, and Xk be any point 
in Bk for 2 < k < m — 1. Then by Lemma [2. Ill since Bi ^ B2 ^' ■ ■ ■ B,n-i -^ Bi 
is a cycle in T„, we have ||/(xfe) — Xfe+i|| < £^ for 1 < k < m. Hence, p is e^-chain 
recurrent. D 

Note this theorem completes the proof of Thcorcm ll.il 

A first modification of this theorem is that if wc consider S' — in the hypothesis, 
we can conclude 3? C 23,i C Ji{e'). We need this observation in Section O to justify 
the process of eliminating "23o-escaping" boxes. 

Next, extending from the theoretical / to the practical F, we immediately obtain: 

Corollary 2.16. Suppose T„ is an (Sn,Sn)-box chain model of F , and F„, e^, (5^, 
and 23„ are as in Theorem \2.1^ 
Let 

Sn{F)^ sup sidelength[F{Bk) - Hull{f{Bk))]. 

BfcGV(r„) 
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Then 

3?(,5;)cS„c3^(e; + s„(i^)). 

We cannot compute Sn{F) exactly, but simply strive to design an implementation 
to minimize it. We shall not discuss that matter in detail in this paper, but refer 
the reader to the interval arithmetic resources listed in Section ^ Away from the 
boundaries of machine precision, it is reasonable to assume s„ (F) decreases as £„ 
decreases, and that s„(F) is small compared to e'„. In fact, for a naive interval 
extension of /, s„(F) = 0(e„). 

We can also immediately apply the theorem to the SCC's. 

Corollary 2.17. Assume the hypothesis of Theorem \2.15[ 

(1) Let Ji'(S') be any 6' -chain transitive component. Then there is a box chain 
transitive component, T' , such that 31' (S') C 23(r'). 

(2) Let V be any box chain transitive component of T. Then there is some 
e'-chain transitive component, 3i'{e'), such that !B(r') C 3?'(e'). 

Again, in Corollary 12.171 with the weaker hypothesis S' = 0, we can form a 
conclusion analogous to (1) by setting 5' = 0, i.e., ^{S') = Jl. 

Note that more than one (5'-chain transitive component can be contained in 
a single box chain transitive component. Also, an e'-chain component may not 
actually contain any chain recurrent points, hence a strongly connected component 
may not contain any chain recurrent points. We explore this further in Sectional 

For polynomial diffcomorphisms of C^ (with d{f) > 1), the Julia set J is con- 
tained in a single chain transitive component 3?' of 3? (see Theorem lA. 1(1 . Hence 
Corollary 12.171 implies : 

Corollary 2.18. Assume the hypothesis of Theorem \2.15\ Then there is a single 
box chain transitive component V such that J C 13 (F'). 

Thus, one of the box chain transitive components, F', contains J . The others 
either contain attracting (or repelling) periodic orbits, or do not intersect D?. We 
can easily identify the component containing J, for it has by far the most vertices. 

Remark 2.19. Note the box chain construction and all of the estimates and results 
of this section can be reformulated to apply to polynomial maps of C of degree d > 1, 
by simply dropping the j/-component, and setting a = in all of the calculations. 
We use these results to study polynomial maps of C in ^21 ■ 

See Section |S1 for examples of box chain recurrent models generated with the 
program Hypatia. 

3. Improving efficiency using Dynamics 

One result of our interest in the complex case is that we are forced to deal 
with increased computational complexity difficulties, since working in C^ is the 
same for a computer as working in R**. In our experience, the main computational 
limitation is memory usage (even using a computer with 4 GB of RAM). Thus 
we keep memory efficiency in mind when tailoring the box chain construction to 
polynomial diffcomorphisms of C^. In this section, we discuss how to insert into 
the box chain construction two improvements, designed to increase efficiency, but 
without completely negating Theorem 12. 151 In designing these algorithms, we take 
advantage of the dynamics of the map. 
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3.1. Step (i'): Selective Subdivision. Note first that a single step (rather than 
inductive) construction consisting of simply subdividing the initial So into a very 
large grid of very small boxes is conceptually much simpler, but very computation- 
ally inefficient. In addition to the obvious work that subdivision saves, the tree 
structure it creates is useful for quickly computing things such as which boxes in- 
tersect a given set (like the image of another box). This allows us to store each 
graph as an array of vertices, with edges as adjacency lists, with no need to arrange 
the vertices in any particular order in the array. Further, the tree structure created 
by subdivision lends itself easily to be improved in the following manner. 

In trying to approximate a dynamically defined set, it is natural to improve upon 
the basic subdivision procedure by replacing step (i) of the box chain construction 
with the more sophisticated selective subdivision procedure. The philosophy here is 
that we can concentrate our resources on the regions which are most "troublesome" 
by allowing boxes of different sizes. We would ideally define grid boxes so that the 
dynamical behavior in each box is varying by at most a small amount. Thus, we 
would like to refine down to a certain reasonable box size, then somehow select only 
a small fraction of the boxes to be subdivided further, leaving the rest unchanged. 

(i') Sort the boxes in V„ into two sets: V^ is boxes to be subdivided, and V^ 
is boxes to remain unchanged. Equally subdivide the boxes in V^: choose 
TO > 1 and place a grid of m^ boxes inside each box of V^, to obtain the 
collection Wfj_|_;^. Set Wn+i = W^_|_;^ U VJJ. Then each box in Wn+i has 
side length at most e„+i ~ £„, and at least En/m. 

Note in this case the sequence of maximum box side lengths {£«} is nonincreas- 
ing, rather than decreasing. In order to acheive the ideal accuracy in the limit, as 
in Theorem 12. 151 one would have to subdivide all of the boxes once every few steps. 

To implement this procedure, we must somehow choose which boxes to subdivide. 
A first goal is to obtain a box chain recurrent model which separates the Julia 
set from the attracting and repelling periodic orbits. With the goal of "refining 
out" the attracting behavior, we implemented in Hypatia an option to subdivide 
only boxes which seem to be in K'^ . We call this sink basin subdivision. For 
example, one way of detecting boxes in sink basins is to choose boxes such that all 
eigenvalues of a few iterates of the derivative matrix Df are small. In some cases, 
sink basin subdivision does help to separate out the sink dynamics more quickly 
(see Example 15.211 . However, there are interesting examples for which this still does 
not allow us to separate the sink from J (see Example 14.1(1 . Thus, it seems that 
a selective subdivision procedure could be an extremely useful step, however it is 
unclear what are the optimal selection criteria for a given map. 

3.2. Step (i.5): Eliminating 23o-escaping boxes. After performing step (i) or 
(i'), i.e., subdividing the desired level n boxes V„ to obtain a new collection of 
level n + 1 boxes, "W„+i, but before performing step (ii) (i.e., computing a graph 
representing the action of / (or F) on the collection of boxes), we institute an 
additional efficiency improving check: eliminating TiQ-escaping boxes. This is a 
very useful technique which removes much of the work of computing the graph in 
step (ii). The idea is to eliminate boxes whose images eventually lie outside of 150. 

Note by examining the proof of Proposition l2.6l that if /™(i?fe) n 23o = for some 
771 > 1, then Bk contains no points oiJi. Thus we can delete vertex k from the graph 
before computing the edge list. For the Hcnon family, we can also take advantage of 
invertibihty. That is, it follows from Lemma ITU that if /'"'{Bk) C C\'Bo, then Bk 
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contains no points of 31. Thus we can check forward and backward images of boxes, 
deleting any boxes with some image leaving 'Bg, and then in step (ii) compute the 
edge list for the reduced graph. 

(i.5) Given a collection of boxes W„+i, eliminate the 'Bo-escaping boxes. Re- 
maining is a subcoUection of boxes y^n+i (which contains 5Vj. Replace 
W„+i with Wj'^+i. 

Each Sn in the constructed sequence contains 31, so we may call it an (e, 0)- 
approximation to 51. Such a sequence satisfies slightly weaker versions of Theo- 
rem 12.151 and Corollarv l2.17l with S' replaced by 0, and 31{S') replaced by Jl. Note 
however we still use the S factor to build the edges in the graph, thus we still satisfy 
Lemma [2.111 and Lemma [2. 131 fthis is important in ^21 and lllp. 

In many instances, this check eliminates half or even three- fourths or more of 
the new boxes. Since the main computational limitation of Hypatia is in memory 
usage, not having to store edges corresponding to boxes which are eliminated with 
this iteration check is a big gain. Specific data for some examples is in Table |21 in 
Section |S1 



4. Separating J from a fixed sink 

In this section, we calculate a theoretical limitation of the box chain construction 
in the case of a Hcnon mapping, /(x, y) = (x^ + c— ay, x), with an attracting fixed 
point p = {z,z). In particular, we establish Theorem II. HI quantifyiing a box size 
which guarantees that the box chain construction produces a box chain recurrent 
model separating the fixed sink from J, as in Definition 1 1.41 i.e., such that the fixed 
sink lies in a separate box chain transitive component from J. This quantification 
is in terms of a,c, and the eigenvalues, Ai ^ A2, of Dpf, and A = max(|Ai| , IA2I). 

First we find a euclidean disk contained in the sink basin, in Proposition 14.81 
Second, we quantify an annular region about the sink which contains only non s- 
chain recurrent points, in Proposition l4.11l We use this along with the estimates of 
Section|2]to derive Theorem II. 61 Finally, we apply our estimates to the 3-1-map. 

Example 4.1 (The 3-1-map). Recall from Example II .51 that the Hcnon mapping 
fa.c with {a,c) = (.3,-1.17) is an interesting example because it appears to have 
two attracting periodic cycles, one of period three and one of period one. Two 
attracting cycles is not a phenomenon which can occur for the quadratic polynomial 
fc{z) ~ z^ + c. Unfortunately, using our program Hypatia implementing the box 
chain construction, we failed to separate the sinks from J before running out of the 
4 GB of RAM available on our computer. 

In an attempt to find a good box chain recurrent set, we first uniformly sub- 
divided all boxes to obtain a (2^)^ grid on Sq = 3\f(0,2.01), with box side length 
2i?/2^ = 0.03, where R = 2.01. Then we used sink basin subdivision (Section OJ, 
which subdivided about half of the boxes. At this point, the smallest boxes had side 
length 2R/2^ = 0.015. The box chain recurrent model F was composed of 944,000 
boxes and 66,500,000 edges. This used approximately 3.2 GB of RAM, thus it 
seemed we could not subdivide significantly farther. 

We also tried uniformly subdividing to obtain a (2^)"* grid on Sq, then invoking 
sink basin subdivision, twice, to get some boxes as small as above. However, this 
did not significantly decrease the amount of memory used. 
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Figure|2]shows the unstable manifold slice of the box chain recurrent set from the 
(2^)^ grid on So- We were initially surprised that we could not achieve separation 
for this map. This motivated the estimates of this section. 

4.1. A dynamically significant norm. To quantify the dynamical notations of 
interest, we need to start with a metric which respects the dynamics. 



Definition 4.2. Note the Jacobian of / is 

D{x,y)f 



2x —a 
1 



Suppose Ai 7^ A2 are the eigenvalues of Dpf, for the fixed sink p = {z,z). Let 
A = max(|Ai| , IA2I). Note that since p is a fixed sink, |A| < 1. Let {vi,V2} be 
the basis of eigenvectors, where we choose Vj = {Xj, 1). Then A = [vi V2] is the 
change of basis matrix, i.e., if {ei, 62} is the standard basis in C^, then Agj = Vj. 
Let ll'llg be the euclidean norm in C^. Define the norm ||-||^ by |ju||^ := ||A^-'^u|| . 

We show below that / is contracting with respect to ||-||^ in an neighborhood 
of p. First; we show this metric is uniformly equivalent to euclidean, and compute 
the constants of equivalence. 

Lemma 4.3. For all u e C^, C||u||^ < j|uj|g < Z3||uj|^, where C,D are positive 
constants given by 



C 



IA1-A2I 



V2+|Ai| + |A2r 



D=^^ 



A2. 



Proof. Let (x,y) be any vector in C^. To calculate both C and D, we use the 
following observation: 

(1) < (|x| - \y\f implies 2 \x\ \v\ < \xf + \yf . 

Recall ||(a;,y)||g. = || A^-^(x, ?/)|| . Since our eigenvectors are Vj = {Xj, 1), 

1 



Ai A2 
1 1 



and A ^ 



Ai — A2 



1 -A2 
-1 Ai 



First we show that we can set C as claimed: 



|Ai - A2r \\ix,y)\\i = |Ai - A2r \\A-\x,y)\\^ = \x - X2y\' + \-x + X,y\' 

< 2 \x\^ + \y\^ (IA2I' + |Ai|') + 2 |a;| \y\ (IA2I + |Ai|), by triangle inequality, 

< 2 \xf + \yf (IA2I' + |Ai|') + {\xf + \yf){\X2\ + |Ai|), by Equationd 
= \x\' (2 + |Ai| + IA2I) + \y\' (IA2I' + |Ai|' + IA2I + |Ai|) , 

< (|a;|' + |2;|')(2+|Ai| + |A2|), since |Ai|,|A2|<l, so IA2I' + |Ai|' < 2. 

Next we need to establish that ||u||_, < D ||uj|^ = D ||A^^u|| , for all u e C^. 
Since A is invertible, this is is equivalent to: ||Au||^ < D j|Au||^ = D \\u\\^, for all 
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u e C^. The following establishes this latter statement with D as claimed: 



|A(a;,i 



|Aia:: + A2?;|^ + \x + y\^ , 



< 
< 

< 

< 



N'dAil' 



y\^ (lAsT + 1) + 2 \x\ \y\ (|Ai| IA2I + 1) (triangle ineq.), 
I' (|Ai|' + |Ai| IA2I + 2) + \y\^ (IA2I' + |Ai| IA2I + 2) (EquationlU, 

x\ +\y\ ) (A + |Ai| IA2I + 2) , since we defined A = max( I All , IA2I), 

I + 2) , since dct(_D/) = a, so |Ai| IA2I — a. 



i2/n(A 



D 



Remark 4.4. Since the eigenvectors arc (Aj,l), the difference |Ai — A2I is the 
determinant of A. This is small if the angle difference between the eigenvectors is 
small. In this case, C captures that the metric is skewed far from euclidean, so only 
a very small euclidean ball can fit inside a cr-ball. Note also that D is large only 
when the eigenvalues are large. Thus D captures the strength of the contraction. 

4.2. Estimating the size of the sink basin. Now that we know how to convert 
between the two norms, we arc ready to take some measurements in the sink basin. 
To do so, we approximate / by Lpf, the linearization of / at p = (z, z). Recall 

Lpf{x, y) = f{z, z) + £»p/(a; - z, y - z) = (z^ + c - 2z(a; - z) - ay, x). 

We first bound the error between / and L^] in the cr-norm. 

Lemma 4.5. Il\\{x-z,y-z)\\^ ^r, then \\f{x,y)-Lpf{x,y)\\^ < r^ (^-^j . 

Proof. Let {x,y) E C^ be such that \\{x — z,y ~ z)||^ = r, for some r > 0. 

It is easy to compute the quadratic error in approximating / with Lpf in the 
euclidean metric, since f{x, y) — Lpf{x, y) = ((x — z)^, 0). 

We then convert to the cr-norm, using Lemma [4.31 twice . to get: 



z\^ < ^\x- z\'^ + \y- z\' 



\iix~z)\0)\\<^\\iix~zr,0)\ 



±\\(:,- z,y- zX< ^\\{x- z,y^ z)t= [^ 



D 



Next, wc show that in the cr-norm, the linearization moves points closer to p by 
a linear contraction. 

Lemma 4.6. If \\{x - z,y - z)\\^ = r, then \\Lpf{x,y) - (z,z)||^ < Xr. 

Proof. Since p is fixed, Lpf(x,y) — (z^z) = Dpf{x — z,y — z). Now to work with 
Dpf, note that since the columns of A are the eigenvectors of Dpf, we have 



A-'DpfA = 



Ai 





A2 



Using this, and the fact that A = max (|Ai| , IA2I), we get 
\\Dpf{x-z,y-z)\\^ = \\A-'Dpf{x~z,y-z)l 

"Ai 



A2 

< X\\A-\x 



A {x — z,y — : 

- z,y- z)\\ = X\\{x 



z^y- z)\\^ = Xr. 



D 
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From the above lemmas and the triangle inequality, we immediately conclude: 
Lemma 4.7. If \\{x - z,y ~ z)\\^ = r, then \\f{x,y)-{z,z)\\^ < Xr + r^ (^) . 

Now we can estimate the euclidean size of the sink basin. 
Proposition 4.8. Let 

IA1-A2I' 



(2+|Ai| + |A2|)(2 + A2 + |a|)- 

Then the euclidean disk centered at p of radius r^ ~ r(l — A) is contained in the 
immediate sink basin of p. 

Proof. Note C^ / D^ = t. We first show that the tj-disk centered at p of radius Sp = 
(1 — A)(C/Z?^), B(j(p, Sp), is contained in the sink basin. For, if |i(x — z,y — z)||^ = 
r < Sp, then by Lemma UTI \\f{x,y) — (z,z)||^ < Ar + r'^{D'^ /C), and by r < Sp 
and definition of Sp we get Xr + r^{D^/C) < r. Thus / maps the disk B(j(p, Sp) into 
itself, and every point in it closer to p in the cr-norm. Thus this disk is contained 
in the immediate sink basin. 

Now if we use Lemma |4. 31 to convert to a euclidean statement, we see that for 



r„ := s 



C we have rp = SpC = (1 - X)iC^/D^) = (1 - A)t, and Be(p,rp) C 
Do-(p, Sp) C {immediate sink basin}. D 

4.3. Separating the box chain transitive components. We now investigate 
the box chain transitive components, i.e., the strongly connected components of F, 
using their relation to the e-chain transitive components as given in Corollarv l2.17l 
First, given a sufficiently small constant ^, we calculate an annular region A(^, 
contained in the immediate sink basin, in which the contraction toward the fixed 
point causes iterates to move toward p by a distance large enough to block ^-chain 
recurrence, with respect to the c-norm. 

Lemma 4.9. LetO < £, < (1 — A)^C/(4D^). Then, in the a-norm, the (,- chain tran- 
sitive component that contains the fixed sink is separated from the ^-chain transitive 
component of any other invariant set by a distance 0/ (1 — X)C/D^. 

Proof. Define ^1^ by ^l^ — {{x, y): r_ < \\{x — z,y — z)\\^ < ?'+} , where p = (z, z) 
is the fixed sink, and 



C 
"±=2D^ 



((1 - A) ± V(i - A)2 - A^oyc) 



We show that if (xq, j/o) S -^J: it is not ^-chain recurrent with respect to || -jl^. Since 
A^ is of CT- width (1 — X)C/D^, this establishes the lemma. 
Note that r± are the roots of the polynomial 

(2) q{r)^iD^/Cy-{l-X)r + t 

Thus ^ < (1 — A)^C/(4D^) is precisely the condition that needs to hold in order for 
the roots r± to be real, and thus positive. Hence, 

(3) q{r) <0, if re ir-,r+). 

Now we show the contraction of / in A^ is strong enough to block ^-chain 
recurrence. Suppose (xq, yo) G A^, with r = \\ {xq — z,yo — z)\\^. Let n £ Z+ and let 
{(xi,yi),...,(a;„_i,y„_i)} be any points such that \\{xj,yj) - /(a-j_i, yj_i)||^ < £,, 



MODELING CHAIN RECURRENCE FOR HENON MAPS 19 

for < j < n — 1. We show that \\{xn,yo) — f{x,i~i,yn-i)\\cr > ■? by first showmg 
inductively that for < j < n — 1, 

(4) \\{xj,yj)~{z,z)\\^<r, so by Lemma EB \\f{xj,yj) - {z, z)\\^ < Xr + r^ — . 

We have (@J) for j = ahcady. Now let < j < n — 1, and suppose we know Q 
for (xj_i,yj_i). Then 

\\{xj,yj)- iz,z)\\^ 

< ll(a^j:%) - .fixj-i,yj-i)\\a + Wfi^j-i^Vj-i) - (2.2;)||^ (triangle ineq.), 

< ^ + Xr + r^{D'^/C) = q{r) +r <r (choice of {xj,yj), Equations [3 and EJ. 

Thus induction verifies I0J. In particular, |0J holds for j = n — 1 and j — 0. But 
then 

\\{xQ,yo)- f{Xn-l,yn-l)\\„ 

> \\{xo,yQ) - {z,z)\\^ - ||/(x„_i,y„„i) - {z,z)\\^ (triangle ineq.) 

> r- Ar-r2(L>VC) ^ ( - lir) >C (Equations H 12 andEJ. 

Hence, {xo,yo) is not ^-cliain recurrent in the cr-norm. D 

A key component of the proof of Theorem 11.61 is the following quantification 
of ?7 > for which the Ty-chain recurrent set is separating, in a sense parallel to 
Definition O 

Definition 4.10. We call the e-chain recurrent set, 31(e); separating if there are 
two chain transitive components, li^ and 'R'^ (of 3^), which lie in different e-chain 
transitive components of 3?(e). In this case we say 3^(e) separates 31-' and 31'"'. 

Further, we call 31(e) fully separating if it separates every pair of chain transitive 
components of 3^. 

Proposition 4.11. Suppose fa^c is a Henon mapping with an attracting fixed point 
p, with Ai 7^ A2 eigenvalues of Dpf, and X = max(|Ai| , IA2I). Let t be as in 
Theorem \1.6\ 

// 77 > satisfies rj < r(l — A)^/4, then ^{r]) is separating. In particular, the 77- 
chain transitive component containing the sink is distinct from the rj-chain transitive 
component of any other invariant set. 

Proof. We have r ~ C^ /D^. We simply convert Lemma lT^ to cuclidean estimates, 
using Lemma E31 Let 7/ = ^C, so that 

^ - ^^ < 4 D2 " 4 ^' 
as claimed in the statement of the proposition. 

Define the set §,, to be simply the set A^, from the proof of Lemma [4.91 Let 
{xo,yo) G §,,, and let 77, € N and {{xi,yi), . . . , (a;„_i, 7/„_i)} be any points such that 
Uxj,yj)- fixj-i,yj^i)\\^ <77, 0<j<n-l. Then \\{xj,yj) - f{x.i^i,yj^i)\\^ < 
T]/C = ^. Thus as in the proof of Lemma [4.91 we have 

||(a;0,yo) - f{Xn-l,yn~l)\\e > C\\{xo,yo) - f{Xn-l,yn-l)\\„ > C(, = T]. 

Thus, [xq, j/o) is not in ^(rj). Thus there exists a connected set §^ in the immediate 
sink basin which lies in C^ \ 3^(77), hence 3^(77) separates the fixed sink from every 
other chain transitive component. D 
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Now we prove Theoreni ll.6l by applying the estimates of Proposition l4. 1 ll to box 
chain recurrent sets, to quantify which box size e guarantees that the components 
of an (e, (5)-box chain recurrent model F separate the fixed sink from every other 
chain transitive component, in the sense of Definition II. 41 

Recall Theorem 11.61 stated that for M > 1 such that 5 < e/M, and k :— 

[l + l/A/ + max{l,(l-A)V^+2||p|| + |a|}], if e < i (-k + ^/k^ _^ t(1 - A)^) , 
then r is separating. 

Proof of Theorem \1.6l Let 23 = 23 (F) be the region in C^ covered by the vertices 
of the {e,5)-hox chain recurrent model F. First note by Proposition 14. Ill if ?; < 
r(l — A)^/4 there is a connected set §,, = 7l^ in the immediate sink basin which lies 
in C^ \ 3^(77). Now if we can calculate a bound on e so that "BnS^ C ^{rj), then we 
would have !B fl S^ = 0. To maximize the bound on e, set 77 = t(1 — A)^/4. 

Now Proposition 14.81 implies Da{p,Sp) is mapped into itself by /, thus there 
are no edges in F from boxes inside ©^(p, Sp) to those outside of ^a{p, Sp). Since 
§,; C Dcr (p,Sp), the box chain transitive component containing the sink must be 
distinct from the box chain transitive component containing J, whenever the box 
size £ is small enough that 23 n §^ C ^{rj). 

To get 23 n 8,, C 3^(77), first recall that Theorem 12 . 1 51 calculates an e' such that 
23 C !K(e'). The theorem specifies that e' = S + e{r + 1), where r is computed in 
Lemma ITTTl as r = e+{2R' + \a\), so e' = 5 + e^ + e{2R+\a\) + e. Note for Henon 
mappings, we do not need to require e < 1, this was applied to higher order terms. 

By examining the proof of Lemma 12.111 and restricting the estimates of that 
proof to apply only to 23 n S,,, we find we can use a slightly better estimate for r. 
Indeed, let i?+ be a bound on the box-norm, ||-||, of a point in S,,. We calculate _R+ 
below. Also recall that by hypothesis we have 6 < e/M. Set 

ly = e/A/ + e2+emax(l,2i?+ + \a\) +e. 

Then the proof of Lemma ITTTl implies 23 n §,, C IR(i'). li v < rj, then 23 n §,, C 
Jl{iy) C 3^(77). So we just need e small enough that v < r]. 

To compute R^ , we use that S,, = A(^ in the cr- norm is centered at the point p, 
and has outer radius 

^+ ^ 2& ((^ " ^) + v/(l-AP-4e^Vc) . 

But at the maximum 77 = ^C, the discriminant is zero, so we have r+ < 252(1 — A). 
Converting to the euclidean norm, we have a bound of Dr+ < ^(1 — A) = \/t{1 — 
A)/2. Now since the box-norm is less than euclidean, we get that if {x,y) £ §,,, 
thcn||(x,7/)|| < V7(l-A)/2 + |b||=:i?+. 

Then ly == s"^ + s[l/M + max{l, (^^(1 - A) + 2 \\p\\) + \a\} + 1], hence notice 
ly = e"^ + en. Now to bound e so that ly < ?/. let q{e) ~ e"^ + en — 77. Then the 
roots of q are (— k ± y^K? + 477)/2. Both roots are real, with one positive and one 
negative. We seek e > small enough that q{e) < 0, which is the same as e smaller 
than the positive root (— k + ^/k^ + ^'q)/2. Since 77 = r(l — A^)/4, this is exactly 
the bound on e claimed in the statement of the theorem. D 

Remark 4.12. In the above proposition, note one does not need all boxes in 23 
to be of size e, but rather just the boxes in the immediate sink basin, computed 
in Proposition l4.8l Thus a selective subdivision procedure targeting the sink basin 
could be advantageous for speeding up separation. 
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Table 1. Constants for sink/ J separation estimates for Example 14. II 



p 


= 


(-0.612,-0.612) 


Ai 


= 


-.885 


A2 


— 


-.34 


A 


~ 


.885 


T 


= 


0.029871571 


t(1 - A) 


= 


0.0034352307 


K 


= 


2.5448759 


V 


= 


9.876288 x 10"^ 


e 


< 


3.880793 X 10"^ 



4.4. One dimension. All of the work of this section applies to Pc{z) — z^ + c in 
the case of a fixed sink p with multiplier A = |P^(p)|. In this case, we do not need 
the cr-norm, so take t = C = D = 1. Then the disk ID)e(p, (1 — A)) is in the sink 
basin and for S < (1 — A)^/4 the set 

As = {z:r^ < \z - p\ <r+}, where r± - i ((1 - A) ± ^(1 - A)^ - 4^) , 

is in C \ 31(6). In order to guarantee separation of J from the sink, for 
K={1 + 1/M + (1 - A) + 2 \p\)) and 77 = (1 - \f/A, 



we need boxes of side length e satisfying e < ( — k + \J ti^ + Arj) /2. 

4.5. The S-l-map. We now apply our estimates to the 3-1-map, to determine how 
small the boxes would need to be to produce a box chain recurrent set separating 
J from the fixed sink. The results of this calculation are shown in Tabled 

Thus, a box side length less than 4 x 10~^ would guarantee that the box chain 
construction separate the fixed sink from J. However, this is several orders of 
magnitude smaller than the best we could compute with current resources (0.015). 
Also, note that the guaranteed euclidcan disk contained in the sink basin is only of 
radius 0.0034. But visually inspecting this example suggests the immediate basin 
is much larger. This suggests a computer should be able to rigorously separate the 
dynamics of J from the sinks, but more sophisticated techniques are needed. 

5. Examples 

In this section we present several examples of applying the box chain construction 
to Hcnon mappings. Every example uses the procedure of eliminating ®o-escaping 
boxes, and some examples use sink basin subdivision. We also present an example 
of the box chain construction applied to a polynomial map of C 

The computations described in this section were run on a Sun Enterprise E3500 
server with 4 GB of RAM and four processors, each 400MIIz UltraSPARC (though 
the multiprocessor was not used). ^ When computations became overwhelming, 
memory usage was the limiting factor. 

To measure the accuracy of an (e,(5)-box chain recurrent model of /, F, we 
compute the bounds e' and 5' given in Theorcm l 1 . II such that 'Jl{6') C 23(r) C 3?(e'). 
As we discussed in Corollary 12.161 the maximum error s{F) between the ideal 



The server was obtained by the CorneU University Mathematics Department through an NSF 
SCREMS grant. 
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r,{z) = z'^ — Sa^z + c, 
with c = —.19 + l.li, a — O.li, and boxes from a 2^° x 2-'^° grid on 
[—2, 2]^. Here IR is J and a 4-cycle. The black regions form the box 
chain transitive component of the 4-cycle. The origin lies in the 
center of the largest black region (in the center of the figure) . The 
box chain transitive component containing J is shaded two tones, 
to heuristically illustrate J. The distance from the origin to the 
end of any of the three limbs of J is approximately 1.3. The small 
shaded islands spread between J and the 4-cycle neighborhood arc 
several box chain transitive components which do not intersect 31, 
but arc contained in some Ji{e'). 

function / and an implemented interval extension of it, F^ evaluated for any box 
in the model, must be added to e' to get the corresponding result for the actually 
computed model. However since at worst s{F) is still much less than e', we neglect 
discussing this factor in the examples presented below. 

Table 121 contains more detailed data for all of the examples of this section. 

5.1. Polynomial maps of C Recall that all of our results can be reformulated for 
polynomial maps of C of degree d > 1. For example, consider the cubic polynomials 
Pa.c{z) = z^ — Sa'^z + c. One can check that R' = 2 suffices for \c\ < 2 and 
|a| < \/2/3. Below we describe an example of an interesting box chain recurrent 
model of a cubic polynomial, computed with Hypatia. 

Example 5.1. The cubic polynomial Pa.c{z) — z^ — 3a^z + c, with c — —.19 -f- 
l.li,a = O.li, has an 31 consisting of J and a 4-cyclc. Figure 13 shows the box 
chain recurrent set which we computed using Hypatia. To get this approximation, 
we uniformly subdivided to obtain boxes from a 2^° x 2^° grid on [—2,2]^. Note 
there are several small box chain transitive components spread between J and the 
4-cycle, which cannot contain any points of 31, but by Corollary 12.171 these box 
chain transitive components are contained in some 3?(e')- 
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5.2. Drawing pictures of box chain recurrent sets for Henon mappings. 

Since applying the box chain construction is an iterative process, we need feedback 
after each step (n) to help us determine when to halt. In this paper, we consider 
a sucessful construction to be one which arrives at a fully separating box chain 
recurrent model (recall Definition II. 4|l . 

To check heuristically whether a box chain recurrent model T is separating, we 
take advantage of the fact that the unstable manifold, W'^{p), of a saddle fixed 
point, p, can be (approximately) explicitly parameterized by the plane (details are 
given in Section |A.3|I . and we sketch the slice of 23 lying in some W'^{p). 

This parameterization identifies p with the origin, and conjugates the Henon 
map to multiplication by the unstable eigenvalue. Hence a parameterized picture 
of J n W"(p) is invariant up to scaling by this eigenvalue, so the restriction to any 
neighborhood of p, [^(p, r), shows a fundamental domain. A box chain recurrent 
set 23 contains boxes of a certain size, so a parameterized 23 n W'^{p) is not quite 
invariant under scaling, but for many choices of r, a sketch of 23 n W"{p) n 'Nip, r) 
will show us whether the box chain transitive components have separated J from 
the sinks, and give an idea of how well 23 is approximating J, especially compared 
to a picture generated the same way with larger boxes (on a previous step). 

To determine the coloration of a pixel, we check whether the pixel intersects 
some boxes in 23. The simplest picture would entail coloring a pixel black if it hits 
23, or white if not. Going one step further, we illustrate the multiple box chain 
transitive components of T by using multiple shades of grey. However, since our 
picture is a parametrization of a one complex dimensional manifold which does not 
line up with the axes in C^, a pixel may hit more than one box, and in more than 
one box chain transitive component. So, we use a color palette which has different 
shades of grey for each box chain transitive component, and a distinctive shade 
(like black) if the pixel hits more than one box chain transitive component. 

We can also enhance our sketch in order to see how well a given box chain 
recurrent set approximates Ji. Theorem lA.il f Section IX!2|) states that J C 01 C K, 
and if W"{p) is the unstable manifold of any saddle periodic point, p, then clearly 
W"{p) C K~ . Hence we slightly lighten the pixels which seem to be in K^; i.e., 
the center point of the pixel stays small after several iterates of the map. In this 
way we can check visually how well 23 approximates J. 

Note that since W^ip) C K~ , the attracting periodic orbits of the map do not 
lie in VF"(p). However, their basins of attraction can intersect W"{p). Thus in 
sketching box chain recurrent sets in this manner, we usually sec part of the box 
chain transitive components corresponding to the sink orbits. 

5.3. Complex Henon mappings. Recall from Section ^ that if fa.c is a Henon 
mapping with a sufficiently small and c is such that the polynomial Pdz) = z^ -\- c 
is hyperbolic (thus stable under perturbation), then /|j is topologically conjugate 
to the function induced by P on the inverse limit lim^( J, P) Ql4|). In this case, we 
say that / is described by P, or simply that / exhibits one dimensional behavior. 

Example 5.2 (The alternate basilica). Recall from Example II .31 that the Henon 
mapping fa^c with c = —1.1875, a = .15, seems to be one of the simplest Henon 
mappings which does not exhibit one dimensional behavior. Computer evidence 
suggests that 31 consists of J and an attracting two-cycle. We found the most ef- 
ficient box chain construction for separating the sink from J was to subdivide 23o 
initially with a (2®)^ grid, then perform sink basin subdivision. In this way we 
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Figure 4. In these pictures, fa.c is the Henon mapping with a = 
— .74, c = —2.75. This map appears to be a horseshoe, thus 51 = J 
appears to be a Cantor set. On the right is a parametrization of 
W"(p), where black is W^{p) f] K~^. On the left is a box chain 
recurrent set with boxes of size 2R' /2^, for R' = 2.84. 

acheived the left side of Figure ^ (shown in Section^. To improve the approxima- 
tion, we performed sink basin subdivision once more, and obtained the right side 
of Figure n (shown in Section^. In these figures, note the box chain transitive 
components skirting the inner edge of the more refined S'. These would not be 
present for smaller box size, i.e., these components do not contain any points of 51, 
though they are contained in some 3^(e'). 

Next we examine a different type of diffeomorphism: a horseshoe. 

Definition 5.3. A horseshoe is a diffeomorphism / such that f\j is topologically 
conjugate to the left shift operator on E2 (the symbol space of bi-infinite sequences 
of O's and I's). The horseshoe locus in the Henon parameter space, !K, is the 
hyperbolic component of parameter space containing the set of horseshoes. 

Since the full left shift on E2 is the inverse limit of the one-sided left shift, 
horseshoes exhibit one dimensional behavior. All horseshoes are topologically the 
same as Smale's horseshoe, i.e., the Julia set is a Cantor set, and the dynamics are 
well- understood . 

However, pictures of the Henon parameter space produced by SaddlcDrop (\l'd\) 
suggest that the topology of the horseshoe locus is quite complicated. There are 
intriguing conjections about the horseshoe locus, which motivate the study of com- 
plex Henon horseshoes. Below we describe one example of applying the box chain 
construction to a horseshoe. 

Example 5.4. The Henon mapping with c = —2.75, a = —.74 appears to be a 
complex horseshoe — though parameters c, a are real, the horseshoe for the map of 
C^ is not contained in R^. For horseshoe diffeomorphisms, there is no sink. Thus 
there is no clear option for selective subdivision. We simply had Hypatia uniformly 
subdivide the boxes, to obtain a box chain recurrent model consisting of boxes from 
a (2^)^ grid on Sq = 3\f(0, 2.84). Since J is a Cantor set, it is easiest to visualize 
J via the FractalAsm picture shown on the left side of Figure 0] The box chain 
recurrent set is shown on the right side of Figure 0] Actually, this cantor set was 
sparse enough that we were able to refine uniformly to a (2^")^ grid, but at this 
point the picture was nearly impossible to see, since the set is very small. 

5.4. Real Henon mappings. The Henon mapping is widely studied as a dif- 
feomorphism of R^, with a,c real parameters. In fact, for some complex Henon 
mappings with real parameters, 51 lies in R^, and the complex dynamics can be 
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Figure 5. The Henon mapping fa,c with {a,c) = (—.25,-3) ap- 
pears to be a real horseshoe, thus 3^ = J C M^ is a Cantor set. 
On the right is a sketch of the stable and unstable manifolds of a 
saddle fixed point. W^{p) U VF"(p), shown in the region [—7,7]^. 
On the left is a box chain recurrent set 23 with boxes of size 2i?'/2^, 
where R' = 2.56. In fact, 23 is a neighborhood of the intersection 

p^«(p)nM^"(p). 



completely described by studying the restriction to M^ (see 0). Observe that all of 
the results of this paper apply immediately to the real setting. We implemented the 
box chain construction for real Henon mappings in Hypatia; below is one example 
of a box chain recurrent set for a real Henon mapping. 

Example 5.5. The Henon mapping fa,c with {a,c) = (—.25,-3) appears to be a 
real horseshoe, in that Jl = J seems to lie in M^. Thus a parameterized unstable 
manifold picture of the Julia set would simply show a Cantor set lying in the real 
axis. With Hypatia we constructed a box chain recurrent set using boxes from a 
(2"^ X 2^) grid on K(0, 2.56) C M^. This box chain recurrent set (shown in R^) is on 
the left in Figure [S] 

The reader familiar with the study of real Henon mappings may recognize that 
the box chain recurrent set (on the left of Figure [SJ appears to show an approxima- 
tion to the intersection of the stable and unstable manifolds, VF"(p) and ^{p) for 
some saddle periodic point p. Indeed, the right side of Figure [S] shows a sketch of 
such manifolds. This observation holds because for any Henon mapping, W^(j)) C 
J- and W'ip) C J+ (see Section EJ. Thus J ^ J+ D J- = W'ip)nW{p). 

Appendix A. Background on Henon mappings 

A.l. The Henon Family. Polynomial diffeomorphisms of C^ necessarily have 
polynomial inverses, thus are often called polynomial automorphisms. Friedland 
and Milnor ([H]) showed that polynomial automorphisms of C^ break down into 
two categories. Elementary automorphisms have simple dynamics, and are polyno- 
mially conjugate to a diffeomorphism of the form {x, y) i-^ {ax+h, cy+p(x)) (p poly- 
nomial, a, c ^ 0). Nonelementary automorphisms are all conjugate to finite compo- 
sitions oi generalized Henon mappings, which are of the form f{x, y) ~ {p{x)—ay, x), 
where p{x) is a monic polynomial of degree d > 1 and a ^ 0. 



Table 2. Data for the box chain models construeted in the Examples of Section [S] Here F' denotes the box chain 
transitive component of F which contains J, and the box grid depth for a box is the number n such that the box 
is of size 2i?'/2". If a box chain model contains boxes of multiple sizes, then multiple box grid depths are listed. 



Recall Example 15.11 is of a cubic polynomial Pa,c{z) 
fa,c{^' y) = (a;^ + c - ay, x). 



3a z + c, and the other examples are of Henon mappings 



<; 

CO 
CO 



Example 


lO 


lOl 


lOI 


E21 


Figure 


n 


m 


HI 


H 


[5] 


params. c 


-.19+ l.li 


-1.1875 


-2.75 


-3 


a 


O.li 


0.15 


-.74 


-.25 


sink period 


4 


2 


N/A 


N/A 


e's.t. Sc3?(£' + s(F)) 


.028 


0.30 


0.30 


0.68 


0.26 


S's.t. Ji{S') c S 


1.5 X 10-*^ 


6 X 10-" 


3 X 10-" 


1.2 X 10-=' 


6.3 X 10-" 


R' 


2.1 


1.9 


2.84 


2.57 


box grid depth, n 


11 


6,7 


6,7,8 


6 


7 


box size 


0.002 


.06, .03 


.06, .03, .015 


0.09 


0.04 


# T boxes original 


98 


184 


682 


60 


20 


(1000s) S-escaping 





116 


417 


43 


15 


T size boxes 


98 


68 


265 


17 


4 


(1000s) edges 


1,300 


2,500 


12,400 


830 


120 


F' size boxes 


60 


53 


182 


10 


2.4 


(1000s) edges 


780 


2,500 


7,800 


500 


75 


runtime (min.) 


<5 




20 


< 1 


< 1 


RAM (MB) 


<200 




900 


20 


20 



MODELING CHAIN RECURRENCE FOR HENON MAPS 27 

To clarify the situation, one can define a dynamical degree of a polynomial au- 
tomorphism of C^. If deg{f) is the maximum of the degrees of the coordinate 
functions, the dynamical degree is 

d = dif) = Km {deginy^-. 

n — 'oo 

This degree is a conjugacy invariant. Elementary automorphisms have dynamical 
degree d = I. A nonclcmentary automorphism is conjugate to some automorphism 
whose polynomial degree is equal to its dynamical degree. Without loss of gen- 
erality, we assume such / are finite compositions of generalized Henon mappings, 
rather than merely conjugate to mappings of this form. 

Thus, the quadratic, complex Henon family fa,c{x,y) = [x^ + c — ay,x) rep- 
resents the dynamical behavior of the simplest class of nonelemcntary polynomial 
automorphisms; those of dynamical degree two. Wc shall state results assuming / 
is a polynomial diffcomorphism of C^ with d{f) > 1, and often concentrate on the 
illustrative example of fa,c{x, y) = {x" + c — ay, x). 

A. 2. Invariant sets of interest. The chain recurrent set. 'R. and the Julia set. J, 
are both attempts at locating the points with dynamically interesting behavior. 31 
can also be decomposed into components which do not interact with one another. 

First we recall the key concepts associated to chain recurrence. An e-chain of 
length n > 1 from p to q is a sequence of points {p = xi, . . . ,Xn = q} such that 
\fixk) — Xk+i\ < e for l<i<n-'l. A point p belongs to the e-chain recurrent 
set, 31(e), of a function / if there is an e-chain from p to p. The chain recurrent set 
is 31 = ne>o3^(e). A point q is in the forward chain limit set of a point p, 3l{p), if 
for all e > 0, for all n > 1, there is an e-chain from p to g of length greater than n. 
Put an equivalence relation on 31 by: p ^ q ii p Cz 3l{q) and q G 3l{p). Equivalence 
classes are called chain transitive components. We can define 31(e) (p) and e-chain 
transitive components analogously. For case of notation, we will sometimes refer 
to 31 as 3i{e) for e = 0, or 31(0). Note that 31 is closed and invariant, and if e < e', 
then 31 c 31(e) c 3l(e'). 

Chain recurrence is quite natural to rigorously study using a computer. A box 
chain recurrent model F is an approximation to the dynamics of / on 31, and the 
connected components of F are approximations to the chain transitive components. 
This is made precise in Section |5] 

For a polynomial map / of C, the filled Julia set, K , is the set of points whose 
orbits are bounded under /; the Julia set, J, is the topological boundary of K. For 
a polynomial diffcomorphism / of C^, there are corresponding Julia sets: K'^{K~) 
is the set of points whose orbits arc bounded under f{f~^) and K ~ K^ n I\^ is 
called the filled Julia set; J* = dK^ (the topological boundary) and J = J^ Ci J~ 
is called the Julia set. The Julia set can be easily sketched by computer, and is also 
extensively used in formulating theoretical results for complex Henon mappings. 

Bedford and Smillie show the following relationships between 31 and J 

Theorem A.l (!$ )• Let f be a polynomial diffcomorphism o/C^, with d{f) > 1. 

(1) Then J <z3l <Z K and J is contained in a single chain transitive component 
of 31. 

(2) Assume further that \detDf\ < 1. Let Oj for j — 1,2,... denote the sink 
orbits of f . 
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(a) Then 'R is the set of bounded orbits (in forward/backward time) not 
in punctured basins, where if p is a sink, the punctured basin of p is 
W^{p) — p, and 

(b) the chain transitive components are the sink orbits, Oj, and the set 

A. 3. Drawing Meaningful pictures for maps of C'^. Filled Julia sets are the 
invariant sets which can be easily sketched by computer, on any two-dimensional 
slice. Hubbard has suggested the following method for drawing a dynamically sig- 
nificant slice of the Julia set of a Hcnon mapping, by parameterizing an unstable 
manifold. This method has been implemented by Karl Papadantonakis into a pro- 
gram called FractalAsm, available for download at J]. We also use this method to 
draw pictures of collections of boxes, see Sections ^ and O for examples. 

Let / be a diffeomorphism of C^. If p is a periodic point of period to, and 
the eigenvalues X, p of Dpf™ satisfy |A| > 1 > |/i| (or vice-versa), then p is a 
saddle periodic point. The large (small) eigenvalue is called the unstable (stable) 
eigenvalue. If p is a saddle periodic point, then the stable manifold of p is W'^{p) = 
{q: d{f'^{q),f'^{p)) ^ as n ^ oo}, and the unstable manifold of p is W'^{p) = 
{q: d{f~'^{q), f~"'{p)) — > as n — > oo}. If p a saddle periodic point of /, then 
W^{p) {W^{p)) is biholomorphically equivalent to C, and on W'^{p) {W'^{p)), f is 
conjugate to multiplication by the unstable (stable) eigenvalue of Dpf . 

Let / be a Henon mapping. When \a\ ^ 1, except on the curve of equation 
4c = (1 + a)^, the map fa.c has at least one saddle fixed point, p, (^Sl)- The 
unstable manifold W'^{p) has a natural paramctrization 7: C ^ VF"(p) given by 

7(z) = lim 7„(z) = lim /™ (p+ t^vi 

where Ai is the unstable eigenvalue of Dpf and vi is the associated eigenvec- 
tor. This paramctrization has the property that f{'y{z)) = 7(Az), and any two 
parametrizations with this property differ by scaling the argument. 

To parameterize VF"(p), we approximate 7 by some g = ^n in a region in the 
plane: B^{z^x + iy: a<x<b,c<y< d}. Observe that since W^{p) C K^, 
to sketch K in VF"(p), we need only sketch /\+. Thus for each pixel Z G B, ii 
f'^{g{Z)) is bounded by some R for all n < N, we guess g{Z) £ K~^ and color Z 
black. Otherwise, color according to which iterate f^{g{Z)) first surpassed R. 

Appendix B. Rigorous Arithmetic 

On a computer, we cannot work with real numbers; instead we work over the 
finite space F of numbers representablc by binary floating point numbers no longer 
than a certain length. For example, since the number 0.1 is not a dyadic rational, 
it has an infinite binary expansion. Thus the computer cannot encode 0.1 exactly. 
Interval arithmetic (I A) provides a method for maintaining rigor in computations, 
and also is natural and efficient for manipulating boxes. The basic objects of I A 
are closed intervals, [a] = [a, a] E IIK, with end points in some fixed field, K. An 
arithmetical operation on two intervals produces a resulting interval which contains 
the real answer. For example, 

[a] + [b] := [a + b,a + b] 

[a] - [b] := [a - 5, a - b] 

Multiplication and division can also be defined in lA. 
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Since an arithmetical operation on two computer numbers in F may not have a 
result in F, in order to implement rigorous lA we must round outward the result of 
any interval arithmetic operation, e.g. for [a], [b] £ IF, 

[a] + [b]:=[ia + bi,]a + b]], 

where |a;| denotes the largest number in F that is strictly less than x (i.e., x 
rounded down), and ta^T denotes the smallest number in F that is strictly greater 
than X (i.e., x rounded up). This is called lA with directed rounding. 

For any a; G R, let Hull(x) be the smallest interval in F which contains x. That is, 
if X e F, then HuU(a;) denotes [x,x]. If a; £ R \ F, then HuU(a;) denotes [i-xi , ta;t]. 
Similarly, for a set 5 C M, we say Hul^^*) for the smallest interval containg S. 
Whether Hull(5') is in IM or IF should be clear from context. 

In higher dimensions, lA operations can be carried out component-wise, on in- 
terval vectors. So if a; € M", then IIull(a;) ~ IIull(a;i) x •• • x IIull(a;„), and if 
S" C R", then Hull(S') is the smallest vector in IF" (or IR") containing S. Note to 
deal with intervals in C" we simply identify C" with R^". Thus a box in C = R^ 
is an interval vector of length two. Our extensive use of boxes is designed to make 
I A calculations natural. 

To compute the image of a point x under a map / using lA, first convert x to 
an interval vector X = IIull(a;), then use a (carefully chosen) combination of the 
basic arithmetical operations to compute an interval vector F(X), such that we are 
guaranteed that Hull(/(a:)) C F{X). 

Definition B.l. Let /: R" -^ R" be continuous. An interval extension of /, 
F = F{f), is a function which maps a box B in R" to a box F{B) containing f{B), 
i.e., F{B) D Hull(/(B)). 

Usually, we would like F{B) to be as close as possible to Hull(/(_B)). We shall not 
discuss how to find the best F. However, in this paper we do distinguish between 
F{B) and IIull(/(_B)) when stating our results, to make clear the difference between 
theory and practice. 

Each time an arithmetical calculation is performed, one must think carefully 
about how to use lA. For example, lA is not distributive. Also, it can easily create 
large error propagation. For example, iterating a polynomial map or diffeomor- 
phism on an interval vector (which is not very close to an attracting period cycle) 
will usually produce a very large interval vector after only a few iterates. That is, 
if i? is a box in R^ = C, and one attempts to compute a box containing fl^{B), for 
friz) ^ z^ + c, by: 

for j from 1 to 10 do 
B = F^B) 

then the box B will likely grow so large that its defining bounds become machine 
oo, i.e., the largest floating point in F. Similarly, if / is a Hcnon mapping, one 
would also never want to try to compute Db^ / o • • • o Db^ f o Dsofiu), for a vector 
u e C^, since the entries would blow up. 

We use lA for all the rigorous computations in the computer program Hypatia. 
The lA routines were all provided by the PROFIL/BIAS package, available at [22| . 
For further background on interval arithmetic, see pi llfilITT] . 
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